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A  Solution  of  th«  8yrtea  of  Nonlinear  Algebra lo 
Balance  Squat Iona  of  Gee -Mixture  by  •Univao-60* 
Digital  Computing  Maohiaa(* 

-  Yugoslavia  - 

/following  la  tho  translation  of  an  artiolo  by 
X.  Paresanovio  and  J.  Patrle  in  Vaslona,  Ho  4,  1969, 
Bolgradoi  pp  120-123. 


Sunnary  of  tha  Whole  Artiolo 

Solution  of  tho  problen  of  finding  tho  balanoo  of  a  gaaooua 
Mixture  ia  rather  frequently  eneounterod  both  by  ohenieal  onginoora 
and  oheniata  dealing  with  rooket  propellant  research.  This  problen 
is  defined  by  a  nonlinear  eysten  of  algebraie  equations  the  computa¬ 
tion  of  whioh  is  a  lengthy  and  tiresome  numerical  job.  The  numerical 
solution  of  such  a  problem  by  desk  computing  machine  requires  several 
days  of  work  for  each  separate  solution  and  is,  of  oourse,  subject 
to  human  errors  made  by  the  calculator  himself.  Zt,  therefore, 
became  apparent  that  the  use  of  digital  oomputing  machine  was 
indispensable.  So  far,  this  problem  has  been  tackled  by  high  memory 
capacity  oomputing  machines,  but  since  many  of  the  engineers  dealing 
with  it  cannot  have  such  me oh In es  made  readily  available,  the  author 
believes  that  they  should  be  given  a  method  of  solving  the  problem 
by  employing  the  optimum  efficiency  of  low- capacity  digital  computing 
machines.  A  complete  solution  is  given  here  for  all  components  of 
the  solution,  this  being  the  most  oomplioated  part  of  the  problem. 

No  numerical  difficulties  should  be  experienced  by  using  the 
j  calculated  compensate  for  plotting  tha  (i,  a)  T-diagrama. 

(Oonoluding  Part)  * 

Tha  carda  were  organised  in  such  a  way  that  it  was  not 
naeasaary  to  transfer  information  from  one  oard  to  the  other,  nor 
waa  it  necessary -to  order  the  oarda.  This  means  that  once  the. cards 
are  ready  and  tha  machine  starts  eoaputlqg,  the  process  continues 
until  the  whole  work  is  finished.  Three  to  four  corrections  were 

AaMjpi.JLTin  Mfrtipu  tot  BragiAB  ,he**I*r»  con*ifyct*d  . 
m  £  This  is  a  oonoluding  part  of  the  artiole.  The  proceeding  part(s) 
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•uch  a  way  that  it  permitted  an  arbitrary  number  of  corrections.  If 
one  performs  k  oorreetiona  to  achieve  oaa  solution,  ona  has  to  taka 
k  ♦  1  oarda.  Note  that  only  the  first  eard  should  carry  (0/90)  and 
the  assuaed  raluas  Is  and  4  . 

The  problem  of  equilibrium  in  the  gaseous  mixture  is  of  suoh 
a  nature  that  ona  has  to  find  components  of  the  composition  of  a 
similar  group  of  problems  share  ai’s  are  fixed  oonstants  while  Ki's 
are  variable  parameters.  Let  us  taka  n  to  represent  the  number  of  j 

similar  problems  and  k  to  denote  the  number  of  oorreetiona  leading  ) 

to  the  solution  of  eabh  problem,  then  the  total  number  of  cards  used  j 
for  the  solution  of  the  similar  group  is  n  x  (k  ♦  1 ).  Only  the  first  j 
oard  among  those  used  has  punched  the  assumed  x?  value.  Sinoe  in 
one  similar  group  there  exist  n  similar  problems  whose  cards  are 
ordered  following  decreasing  or  inoreasing  values  of  Ki,  it  seems 
convenient  to  use  the  corrected,  value  for  xs  from  the  previous 
problem  as  the  assumed  value  for  the  next  problem  of  the  same  group. 
This  means  that  the  oomputed  value  for  xe,  i.e.  x,  obtained  in  the 
first  correction  serves  as  the  assumed  value  for  the  next  correction. 

Assuming  the  value  xff  instead  of  x,  one  introduoes  one  of  the  , 
big  simplifications  into  the  solution  of  this  problem.  In  this  way,  j 
the  solution  is  obtained  much  faster,  the  setting-up  of  the  program 
is  tax  oh  simpler,  and  the  use  of  the  memory  is  much  more  economical. 

One  oan  conclude  from  the  large  number  of  problems  solved  using  the 
UNIVAC-60  that  each  solution  requires  one  minute  of  time. 

Pig.  1  shows  the  flow  diagram  for  program  I.  One  should 
mention  that  during  the  calculation  of  the  first  correction  of  x£ 
one  does  not  erase  the  storages  S2  and  83  which  permits,  as  mentioned 
earlier,  that  the  computed  value  for  xs  from  the  previous  problem 
can  be  used  as  the  assumed  value  for  the  next  problem  belonging  to 
the  same  similar  group  of  problems.  When  computing  the  second, 
third,  etc. 1  correction  one  should  erase  the  storages  82  and  83. 

The  program  selector  PSi  permits  a  selection  of  data  in 
program  steps  5-15  and  together  with  PS2  and  PS4  makes  the  iterative 
process  of  xs  correction  possible,  while  the  PS3  determines  if  one 
has  to  inereaee  or  decrease  the  Revalues  during  the  correction.  On 
Pig.  2  are  indicated  the  connections  between  various  selectors  used 
in  program  I. 


Elements  and  storages  uaad  for  program  I  contain  tha  following 


valuoa. 


[nj]  — -  X  a  assumed 

(N»|  —  2  K, 

IN,]--  0 

|N,  J  —  K, 

(N,,)-— K„  i.e.*  a 

{N„l— a. 

(S,|-—  Z 

|N,|=-  «, 

IN,!— K„ 

IS.I  X2  calculated 

|N,1  —  K,.  i.e. «  1°  3 

jS,|  -  X  Ticalculated 

|N*|  a. 

|S,J=—  y 

Values  taken  off  the  cards  are 

xTass  K-.  K,.  Kw  i(lji  (2) 

while  the  values  from  the  program  plates  are 

aj ;  a^;  a^:  aj  o;  0.001  and  0.001  (5) 


The  flow  diagram  for  program  II  ie  presented  in  Fig.  5* 

The  first  step  in  this  program  should  allow  the  switching  of 
selector  2  if  there  is  (0/90).  If  this  step  is  missing,  it  oould 
happen  that  the  selectors  do  not  ohange  to  the  "selector"  position  j 

_  i 


end  the  program  could  proceed  in  the  wrong  way.  The  result  of  the 
15th  program  step  operation  is  always  negative  and  is  utilized  for 
punching  the  (1/90).  The  program  will  be  executed  beyond  this  step 
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only  If  the  value  for  d  obtained  satisfies  the  required  accuracy. 
Otherwise,  the  machine  computed (4  /p )t  and  the  card  ia  sorted  out* 
Por  the  following  cards,  from  the  group  of  cards  for  a  single 
solution,  the  program  passes  or  does  not  pass  through  the  20th 
program  step  depending  on  the  fact  was  the  (1/90)  punched  on  the 
previous  eard,i.ev  does  the  solution  satisfy  the  required  accuracy 
or  not.  If  the  previous  solution  requires  further  correction,  the 
following  card  is  punched  with  (A/p)t  eaXc.»  d  calc.»  Rnd  (°9o) 

since  the  result  of  the  arithmetic  operation  in  the  20th  program 
step  is  always  negative .  The  connection  with  SSL.  HOLD  is  broken 
when  P Sj  enters  the  operation.  Fig.  4  shows  the  connection  of  the 
selector  11a  and  4  which  are  used  to  maintain  selectors  11b  and  lie 
on  the  selector  side  after  the  TRIP  command. 
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Pig.  4. 


Elements  and  storages  used  in  program  II  contain  the 
following  values 


IN,!  -- 


1 


|N,1=^  K,  1 

iNt]  —  y 

{N.1  «: 

l 

|N4j x*  calculated 

H— (if 

|N,1  K,  <i.e  .0  lap— 5,1=0.004 

[N„|  =►  a.  oi*®»3  * 


IN,. I  -  -  K, 

|N»!  2 

|Ntt|=»  P 
[S»l*=*'  « 


CM 


, 


i 

i 


i 

j 


i 


This  procram  us  as  elements  xe,  Y,  K} ,  Kjt  K^,  S  calc.*  and 
from  the  card,  and  a-j,  aj,  a^,  a,  p,  2,  1,  and  0.004  from  the 
program  plates. 

Fig.  5  indicates  the  connections  between  individual 
selectors  used  in  the  program  II. 

The  flow  diagram  for  program  II  i9  shown  on  Fig.  6.  This 
program  is  utilized  after  the  solutions  are  found  for  x,  y,  and  z 
from  all  problems.  These  solutions  are  punched  on  cards  with  the 
(1/90)  label  and  after  each  such  card  one  should  insert  two  empty 
cards  •which  will  be  punched  with  the  Xi  (i*»1,2,..».»  10) 
solutions.  During  the  first  passage  will  be  calculated  and  punched 
the  values  for  xj,  *4,  X5,  x£,  and  xg,  while  the  second  passage 
will  yield  x1,  X2,  *7,  X9,  xio»  The  program  for  finding  these  ten 
components  is  constructed  after  Equations  (17-26).  Elements  and 
storages  used  contain  the  following  values? 
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[n'HD! 

|N,}  —  2K, 
IN,}  —  y 
IN*)  —  Z 
IN,}  —  Kl0 
IN,}—  K? 
IN,}—  K, 

i 

|N,1—  X2 
|N,)  —  K, 


IN,,)  —  a, 

(Nn|  —  2 
IN,,)  —  a, 

|S,]  =—  X,.  <1.6.3  X, 

fS,)  —  x&.  ri.e.,  xw  (5) 
is*)  —  X*.  <!.«.,  X. 

|S,1—  X,.  I  i.®«  »  X, 

IS,)— X*.  .i.e.o  X, 


Pig.  5 


ebi^isj 


?'r?l 


I 

4(pfl*l|3^62 

*lBB*Bg=»Bia 

jfgg «;  oflt  =»gg 

mu^=j>sa 

olBO:'  &d=i>Bil 


&a-ea=>B0 


X' 

scrf\  3/bo 


^sa»&<l=>&a 

&a|  ^(STI 
TZ± 

*l£fc)«  fea  =>&*} 

«(lia* 

<y{fy<D» 


J  Smr 


Ttm  the  sards  are  utilised  the  values  of  »,  y,  t,  *i,  2X5* 
I7,  I©,  Kio*  *ad  (d/p)*,  aMtl#  the  plates  supplied  a*,  a«,  and  2* 

'  Pig.  7  indicates  the  dense et  ion  between  individual  selectors 

used  in  the  prograa  III* 


Pig.  7. 


VI  -  Conclusion  , 

This  article  presents  the  possibility  of  solving  the  problem  j 
of  equilibrium  in  a  gaseous  mixture  using  digital  computers  with  j 

small  capacities  of  memories.  The  working  equations  are  given  in  a 
convenient  form  which  allows  a  substantial  saving  of  program  steps 
independently  from  the  type  of  digital  computer  used  for  the  solution  i 
of  this  problem.  The  assumption  of  x*  instead  of  x  was  cf  great 
help  during  the  solution  of  working  equation.  After  finding  all  the 
components  of  the  composition  of  the  gaseous  mixture  one  can  consider 
as  solved  the  bulk  of  the  numerical  calculation  required  for  the 
practical  application  of  this  problem.  Further  numerical  work  does 
not  present  any  difficulty* 
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